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Frequency ratio method for seismic modeling of y Doradus stars 
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Abstract. A method for obtaining asteroseismological information of a y Doradus oscillating star showing at least three pul- 
sation frequencies is presented. This method is based on a first-order asymptotic g-mode expression, in agreement with the 
internal structure of y Doradus stars. The information obtained is twofold: 1) a possible identification of the radial order n and 
degree t of observed frequencies (assuming that these have the same £), and 2) an estimate of the integral of the buoyancy 
frequency (Brunt-Vaisala) weighted over the stellar radius along the radiative zone. The accuracy of the method as well as 
its theoretical consistency are also discussed for a typical y Doradus stellar model. Finally, the frequency ratios method has 
been tested with observed frequencies of the y Doradus star HD 12901. The number of representative models verifying the 
complete set of constraints (the location in the HR diagram, the Brunt-Vaisala frequency integral, the observed metallicity and 
frequencies and a reliable identification of n and I) is drastically reduced to six. 
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1. Introduction 

The asymptotic approximation for the solution of linear, isen- 
tropic g-mode oscillat i ons of high radial order w a s firstl y at- 
tempt ed by ISmeversI (I1968I) . iTassoul & T assou 3 ill 968b and 
IZahnl Jl970l) . These initial approaches were developed for the 
special case of a star with a convective core and a radiative en- 
velope. From these first developments until recently, a ll asymp- 
totic s olutions have been carried out by using the ICow ling 
Jl94lb approximation, where the perturbation of the gravita- 
tio nal poten t ial is neglected. 
iTassou 

linn 

investigated the asymptotic representation 
of high-frequency p-modes and low-frequency g-modes asso- 
ciated wit h low-degree spher i cal ha rmo nics. There, the pro- 
cedure of llweins & Smeversl dl968b and IVandakurovl (1 1967b 
was adopted to avoid the difficulties of the so-called mobile 
singularities. This procedure solves different equations in dif- 
ferent regions of the star and gives asymptotic expressions 
for the p and g-modes fo r different stellar structures. Later, 
ISmevers & Tassoul Jl987b used pulsating second-order differ- 
ential equations, expressed in a single dependent variable, to 
make the construction of successive asymptotic approxima- 
tio ns more t r anspar ent. 

Tassou 3 J 1990b developed, for the first time, an asymp- 
totic expression for high-frequency non-radial p modes with- 
out neglecting the Eulerian perturbation of the gravitational 
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potential. This study wa s carri ed out for low-frequency g- 
modes by ISmevers et al. I d 1995b by using a different asymp- 
totic approach. They applied a procedure described by 
iKevorkian & Cole] Jl98lb to a fully radiative star, which was 
extended to a convective core plus a radiative envelope by 
Willems. Van Hool st & Smeversl (11997b . Since the work of 
Tassoul (1980), subsequent studies have improved the second 
order and the eigenfunction descriptions, but the first order 
asymptotic expressions for the freq uency in differe nt physical 
situations was correctly obtained by Tassouj Jl980b . 

The first attempt to obtain information on a real star by 
usin g these g-mode asymp t otic d evelopments was carried out 
by IProvost & B erthomi eul dl986b ada pting the second order 
asymptotic theory of lTassoull (1980) to the special case of 
the Sun. They argued that the study of the g modes of a 
star would provide information about th e physica l condi tions 
of the stellar core. At th e same time, iKawaleJ dl987b and 
iKawaler & Bradlevl \ 1994b started to investigate the internal 
structure of the PG 1 159 hot white dwarf stars through the same 
asymptotic g-mode theory. By using the fully radiative star ex- 
pression, they obtained stellar properties comparing mode dif- 
ferences of consecutive orders. 

Re cently, y Doradus stars have been defined bv lKave et alJ 
(1999) as a class of stars pulsating in high-n, \o\n-{ g modes 
with very low photometric amplitude, lying on or just to the 
right of the cold border at the lower part of the Cepheid insta- 
bility strip. During recent years, tens of candid ate y Doradus 

stars have been observed an d cataloged (Mathias et al. 2004; 

L II 

Martin. Bossi & Zerbi 2003, and references therein). Also, sev- 
eral theoretical studies have been carried out to understand their 
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instability mechanism, since the ^--mechanis m does not ex plain 
the observed modes (see iGuzik et"afll2000t iGrigahcene et alJ 
l2004lDupret et all2004l) . 

Modelling these and other pulsating stars is mainly de- 
pendent on the fundamental stellar parameters such as metal- 
licity, effective temperature, luminosity and surface grav- 
ity as deduced from observations. In the case of the 
well-known five minute solar oscillations, helioseismol- 
ogy has provided important improvements to our knowl- 
edge of the internal structure and evolution of the Sun 
than ks to the asymptotic theory applied to solar p modes 
( Chr istensen-Dalsgaard. Thompson & Goughl 19891) . useful for 
constraining the mass and the evolutionary stage of solar-like 
pulsating stars by using the (Aero, -Do) diagrams. Acr repre- 
sents the average separation between modes of the same de- 
gree and adjacent radial order and Do is related to the small 
separation between <x„ / and o- n -\,ui- This increases the set of 
parameters CP provided by classical observables with the seis- 
mic information 

CP = y(Z,g,L,T eff ;Ac7-o,Do). 
where Z represent the relative metal abundance, g the gravity, 
L the luminosity and T e ff the effective temperature of the star. 

An analogous method for low frequency asymptotic g- 
mode pulsators, such as y Doradus stars, will be given here. 
Additional constraints can be obtained by using information 
inferred from ratios of the observed frequencies. Similarly to 
acoustic modes in the asymptotic regime, our method increases 
the set of parameters CP provided by classical observables by 

CP = CP(Z,g,L,T eff ;^,CI). 

Jo,j 

where / 0; , and f a j represent the observed frequencies and 3 
is the integral of the Brunt-Vaisala frequency defined below. 
Particularly, this procedure will provide an estimate of the ra- 
dial and spherical orders (n, €) of at least three observed fre- 
quencies of a y Doradus star. In addition, the internal stellar 
structure will also be constrained through the knowledge of the 
integral along the stellar radius of the Brunt-Vaisala frequency. 

The paper is organized as follows: In Sec. 2, the first or- 
der analytical asymptotic equation to be used is presented, as 
well as the physical reasons for the choice of this expression. 
In Sec. 3, the procedure and a test of its theoretical consistency 
is discussed. An application of the method to observations of 
the y Doradus star HD 12901 is given in Sec. 4. 

2. Analytical equation for asymptotic g modes 

y Doradus stars are found to pulsate in the low frequency g- 
mode asymptotic regime. There are several analytical forms 
for obtaining the numerical value of the pulsational period as 
a function of the radial order n, the spherical order I, different 
constants and equilibrium quantities, all of them having slight 
differences and are applied to different internal stellar struc- 
tures. 

These stars present a convective core, a radiative enve- 
lope and a small outer convective zone close to the photo- 
sphere. The distribution of convective and radiative zones in- 
side the star produces different propagation cavities, which is 




Fig. 1. Propagation diagram for the theoretical model presented 
in Table \l\ Angular frequencies are given in nondimentional 
units. The horizontal line represent a standard y Doradus fre- 
quency, around 20 pHz (« 2 = 0.0018). 



the key for solving the problem of interpreting the behavior of 
their g-modes in the asymptotic regime. The so-called turning 
points represent the limits of such cavities. They are defined, 
depending on the a uthors, by the zeros of the Brunt-Vaisala 
frequency jTassouli ri980) or by the points of intersection of 
the propagating mode with the Brunt-Vaisala an d Lamb fre- 
quencies, whichever comes first Jshibahashi ll979l) (see Fig.Q. 
Both quantities represent the characteristi c buoyanc y (Brunt- 
Vaisala) and acoustic (Lamb) frequencies JUnno et al.ll989h . 

For the y Doradus internal distribution of convective and 
radiative zones, the better adapted analytical solution is 



1/2 



(«+ l/2)7r" 

given bv lTassoull(ll980l) and obtained under the assumption 
of adiabaticity and non-rotation. 3 is the integral of the Brunt- 
Vaisala frequency, given by 



3 



(2) 



where r a and r/, represent the inner and outer turning points 
respectively. 

Considering low g-mode frequencies in the asymptotic 
regime (typical of y Doradus stars), the inner turning point 
(r a ) is clearly defined by the zero of A^ 2 (Fig.Q. However, for 
the outer layers, small changes in the frequency produce dif- 
ferent turning points (fb) depending on whether the frequency 
intersects first A^ 2 or L\ (see Fig. [3 compared to the zero of 
Af 2 . The problem of considering one of these possibilities re- 
mains unsolved. This theoretical question beyond the scope of 
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20 25 
Frequency (llHz) 

Fig. 2. Relative errors of asymptotic frequencies compared to 
theoretical predictions (e r<cr ) by using Eq^for a complete the- 
oretical oscillation spectrum. See text for details. 

the present work, however an estimate of the error committed 
when considering Eq.^can be calculated. To do so, theoretical 
oscillation frequencies <x„ £ computed from a given model (see 
Table are compared with those obtained from Eq. \l\ Such 
errors, defined as 



|°"».f-<f| 



<r n ,t 



(3) 



are shown in Fig.[2]for 6= 1 and I — 2 modes. As expected, e v 
increases with the mode frequency and hence decreases with 
n, since they depart from the asymptotic regime. In the range 
of frequencies in which y Doradus stars pulsate, these errors 
remain around 1-2%. 

Equivalently, the relative error committed when calculating 
3 values from Eq.^ can be investigated through 



e r ,3 



1 3th - 3" | 



(4) 



where 3? is calculated for each mode provided the frequency 
and the radial and spherical orders of the modes, and 3th repre- 
sents the integral 3 calculated for the theoretical model. Errors 
obtained in this case are of the same order as e r ,<r, and the aver- 
age value of 1% observed in Fig.[2]will be used in Sec. 4. They 
remain lower than typical errors given by the resolution of ob- 
served frequencies and hence, of their ratios. This represents an 
important advantage for the method since no significant addi- 
tional uncertainties need to be considered besides those coming 
from observations. 

Therefore, consider ing all p r eviou s arguments, the analytic 
form (Eq. [0 given by Tassoul ( 1980) will be adopted in this 
work. 

3. Frequency Ratio method 

3.1. Method bases 

In the case of asymptotic acoustic modes, physical information 
can be inferred from frequency differences of adjacent radial 



order modes. This can be done since the model dependence ap- 
pears in the equation as an adding term. As can be seen from 
the analytic form for asymptotic g modes (Eq.Q, it is not possi- 
ble to retrieve physical information from frequency differences. 
However, the role played by such frequency differences can be 
replaced by frequency ratios. 

Let us consider two modes cr ai and cr a2 , with a\ = 
and a2 = («2> (2) respectively. For sake the of simplicity, the 
same mode degree t\ — I2 — I of a non-rotating pulsating star, 
showing adiabatic g modes in asymptotic regime, is considered. 
Their eigenfrequencies can be approximated by Eq.^ which is 
model-dependent through the integral given by Eq.|2] As shown 
in Fig.[0 for typical y Doradus frequency ranges, the possible 
slight differences of the outer turning point location can be con- 
sidered negligible as far as the calculation of 3 is concerned. It 
can thus be approximated as constant (3 a2 = 3 ai = 3), and 
therefore, the ratio of cr ai and cr 0l can be approximated by 



cr ai ^ »2 + 1/2 

<T a2 Hi +1/2 



(5) 



An estimate of the radial order ratios of these frequencies can 
thus be obtained. A value of the integral can also be deduced 
from observations 3,,/™, provided that an estimate of the mode 
degree t is assumed. Constraints on models will come from the 
consistency of: 1) the radial order identification corresponding 
to the observed ratios; 2) their corresponding order i and 3) the 
observed Brunt-Vaisala integral. 

On the other hand, the assumption of equal t for all the ob- 
served modes is not imperative for this procedure. Additional 
information on £ provided by spectroscopy or multicolor pho- 
tometry can b e exploited through the following expression (see 
Tasso ul 19801) 



<r ni A „ "2 + 1/2 Vgjgj + 1) 
cr„ 2 ,f 2 ~ Mi + 1/2 y/i 2 d2 + 1) 



(6) 



which is the extended form of Eq.|5] 

The efficiency and utility of this method depend on the 
number of observed frequencies N/. The larger N/, the lower 
the number of integers verifying Eq. [5] The method becomes 
useful for N/ > 3. 

3.2. Theoretical test 

To test the self-consistency of this technique, consider a 
theoretical numerical simulation of a given star. Assuming 
three theoretical frequencies computed for a 1.5M© typical 
y Doradus model on the Main Sequence, the ratio method is 
applied to obtain an estimate of the radial order n and the spher- 
ical order I of these frequencies as well as their corresponding 
3. If this procedure is self-coherent, at least one of all possible 
solutions will correspond to the selected theoretical frequencies 
and the corresponding 3 value. The theoretical eigenfrequen- 
cies (computed by the code FILOU, see Sec. 4.2) of a typical 
y Doradus model, as well as its 3^ value in fiHz are given in 
Table[T] 

A first test of the ratio method is carried out for three close 
I = 1 modes, for instance (n,cf) = (25, 14.706), (23, 15.922) 
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Table 1. £, n and frequency cr in fj,Hz for a complete spectrum 
of the theoretical model displayed at the top of the table. 



log T cfI log g log ± 



V M 



3.845 


4.045 


0.873 


0.29 


1.5 


839 


I 


ii 


cr(nHz) 


I 


n 


trQxHz) 




30 


12.3052 


j 


22 


16.6262 




29 


12.7176 


j 


21 


17.4038 


j 


28 


13.1589 


j 


20 


18.2397 




27 


13.6375 


j 


19 


19.1344 


! 


26 


14.1553 


j 


18 


20.1369 




25 


14.7063 


1 


17 


21.2964 


! 


24 


15.2893 


j 


16 


22.6059 


j 


23 


15.9219 


j 


15 


24.0274 


2 


53 


12.1277 


2 


40 


16.0399 


2 


52 


12.3621 


2 


39 


16.4502 


2 


51 


12.5999 


2 


38 


16.8789 


2 


50 


12.8444 


2 


37 


17.3283 


2 


49 


13.1012 


2 


36 


17.8010 


2 


48 


13.3737 


2 


35 


18.2982 


2 


47 


13.6612 


2 


34 


18.8212 


2 


46 


13.9616 


2 


33 


19.3752 


2 


45 


14.2724 


2 


32 


19.9675 


2 


44 


14.5931 


2 


31 


20.6010 


2 


43 


14.9264 


2 


30 


21.2734 


2 


42 


15.2772 


2 


29 


21.9839 


2 


41 


15.6484 


2 


28 


22.7434 



(7) 



and (20, 18.240). The corresponding frequency ratios are 

CTi (Ti (T\ 

— = 0.9236, — = 0.8729, — = 0.8063 

cr 2 cr 3 cr 3 

which constitute our assumed observed frequency ratios. As in 
this theoretical test the radial order of selected frequencies is 
known, it is possible to study the accuracy of Eq.|5]by compar- 
ing the corresponding n ratios: 

m + 1/2 



m + 1/2 



= /Si = 0.9216, p 2 = 0.8723, #, = 0.8039 (8) 



yielding an error of 2x 10" 3 , 6x 10~ 4 and 24x 10~ 4 respectively. 
In Fig. [3] the accuracy of Eq.|5J defined as 



a-, Hj + 1 12 



(9) 



P CTj n t + 1/2 

is shown for all the frequencies of Tabled The maximum de- 
parture of the n ratio from the real frequency ratios is 5 ■ 10~ 3 , 
allowing us to consider this value as the error when applying 
this method. 

3.2.1 . Radial order identification 

If the values of the radial orders are assumed unknown and their 
errors are supposed to be given by the maximum of ep (Eq.[9}, 
i.e. 5 ■ 10~ 3 , ratios of integer numbers are searched for within 
the following ranges 



cr 2 



= [0.929,0.919] 



0.005 
0.004 
0.003 
0.002 
0.001 



mfi* 
0.001 



-0.002 



* * + 



L=l 

L=2 



8 iS 



mi 



12 14 16 18 20 22 24 26 28 
Frequency (liHz) 



30 



Fig. 3. Error ep between the frequency ratios and the corre- 
sponding asymptotic predictions, for the complete spectrum of 
the theoretical model. 

Table 2. Radial number n for each selected frequency, I and 
Brunt-Vaisala integral (in juHz) estimated by the frequency ra- 
tio method for the assumed observed frequencies (20, 18.2397), 
(23, 15.9219) and (25, 14.7063) with £ = 1. The real values of 
the radial and spherical orders are displayed in boldface. 



"1 


"2 


n 3 


e 






20 


23 


25 


i 


831 


0.0017 


33 


38 


41 


i 


1360 


-0.0007 


33 


38 


41 


2 


784 


-0.0007 


41 


47 


51 


2 


971 


0.0003 


21 


24 


26 


1 


702 


-0.0035 


42 


48 


52 


2 


801 


-0.0023 


34 


39 


42 


1 


1127 


-0.0039 


34 


39 


42 


2 


651 


-0.0039 


40 


46 


50 


2 


764 


0.0030 



— = [0.878,0.868] 

0~3 



(10) 



11 

0"3 



[0.811,0.801] 



Integer numbers contemplated here (from 1 to 60) cover 
the range of typical g-mode radial orders shown by 
y Doradus stars. Within this range, there exist several pos- 
sible integer numbers yielding the same value for each ratio 
of Eq. However, the key point consists of searching for 
(«i,«2, «3) sets verifying 



«2 + 0.5 cr; n 3 + 0.5 <t 2 «3 + 0.5 
«i+0.5 0"2 «2 + 0.5 0"3 n\ + 0.5 



^ (ID 
cr 3 



and limiting thus the final number of valid sets. 

Once all possible sets have been obtained, an estimated 
spherical order £ to which ti\, ri2 and n 3 are associated can also 
be identified. To do so, selected sets are compared with theo- 
retical oscillation spectra computed from representative models 
of y Doradus stars. Finally, for each solution (a,n,f) found, an 
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Table 3. Radial number n for each selected frequency, I 
and Brunt-Vaisala integral estimated by the frequency ratio 
method for the assumed observed frequencies (16,22.6059), 
(23,15.9219) and (28,13.1589) with t = 1 (set 1), 
(45, 14.2724), (47, 13.6612) and (50, 12.8444) with I = 2 
(set 2) and (30,21.2734), (45, 14.2724) and (53, 12.1277) with 
1-2 (set 3), See text for details. The real radial and spherical 
orders are displayed in boldface. 



set 


»i 


m_ 


"3 


e 






1 


26 


37 


45 


2 


768 


-0.0001 




33 


47 


57 


2 


971 


-0.0004 




16 


23 


28 


1 


830 


0.0024 




23 


33 


40 


1 


1182 


0.0013 




23 


33 


40 


2 


682 


0.0013 




30 


43 


52 


2 


886 


0.0007 


set 


n\ 


"2 


ns 


e 






2 


46 


48 


51 


2 


850 


-0.0020 




45 


47 


50 


2 


832 


-0.0007 




43 


45 


48 


2 


798 


0.0021 




44 


46 


49 


2 


815 


0.0006 




47 


49 


52 


2 


869 


-0.0033 


set 


»i 


"2 


"3 


e 






3 


22 


33 


39 


i 


1063 


0.0005 




22 


33 


39 


2 


614 


0.0005 




26 


39 


46 


1 


1252 


0.0001 




26 


39 


46 


2 


723 


0.0001 




30 


45 


53 


2 


833 


-0.0004 




34 


51 


60 


2 


941 


-0.0002 




33 


49 


58 


2 


910 


-0.0016 




18 


27 


32 


1 


876 


0.0008 



estimate of 3th is approximately constant within each set) can 
also be deduced through Eq.^ 

In this exercise seven possible sets have been obtained. The 
real n and I (our simulated observations) have been retrieved 
(boldface in Table [3 together with other eight possible iden- 
tifications. In addition, an accurate prediction of the Brunt- 
Vaisala frequency is also obtained. To investigate other distri- 
butions of observed frequencies (always within the asymptotic 
regime), the following cases are studied: 1) modes with very 
different radial order, for instance (n\ - 16, «2 = 23, «3 = 28) 
and 6 = 1; 2) close consecutive radial orders, for instance 
(«i = 45, «2 = 47, «3 = 50) and 1-2; and 3) very different 
radial orders (n\ - 30, «2 = 45, «3 = 53) and 1-2. Following 
the same procedure, sets of 6, 5 and 8 possible solutions are 
obtained. For all of them, the real values of n and i (boldface 
in Table [3), and an estimate of 3 very close to the values from 
the model are obtained. 

These results are very important since they show the self- 
consistency of the method which, in turn, reduces drastically 
the number of representative models of a given y Doradus star. 
Moreover, they illustrate the good accuracy of Eq. ^ f° r this 
kind of study. 



Table 4. Photometric observed frequencies in both cycles per 
day and pHz for the y Doradus star HD 12901, taken from 
AertKetalJiaXMl. 





V 


V 




(c/d) 


CuHz) 


fl 


1.216 


14.069 


fn 


1.396 


16.157 


fm 


2.186 


25.305 



Table 5. Physical parameters of the stars HD 12901 taken from 
the literature. The values of the rotational velocities are also 
listed. 



HD 


T e ir 


log g 


[Fe/H] 


vsim 


Ref. 




(K) 


(dex) 


(dex) 


(km s" 1 ) 




12901 


6996 


4.04 


-0.37 




a 




7079 


4.47 


-0.40 




b 










53 


c 










66 


d 



a: Temp LofjG 1 Rogers 1995) 
fr: |DupretN2002l) 
c: Aerts et al. 120041) 
rf: llvIathiasetalJ<2004 



4. An application: The y Doradus star HD 12901 

Once the frequency ratio method has been shown to be self- 
consistent, a test with real data must be carried out for the 
y Doradus star HD 12901. This star pulsates with three fre- 
quencies (Table 0} in the asymptotic regime (see Sec. [3] for 
more details). 

4.1. Fundamental parameters 

The fundamental parameters of H P 12901 have b een de- 
termined by applying TempLogG (Rogers 1995) to the 
Stromgren-Crawford photometry listed i n the Hauck- 
Mermilliod catalogue jHauck & Mermillio J Il998h . In this 
catalogue, no (3 index value is given f or HD 12901, which was 
instead obtained from Handler ( 1999). The code classifies this 
object as a main sequence star in the spectral region F0-G2. 
The resulting physical parameters are summarized in Table[5] 
Dupret (2002) gives also physical parameters for HD 12901 
based on the Geneva photometry o flAer ts et alJ d2004l) and the 
calibrations given by Kunzli et al ] dl997l) for the Geneva pho- 
tometry of B to G stars. In his study, theoretical stellar models 
for this star do not take into account their surface gravity deter- 
minations, as they are too high, corresponding to models below 
fheZAMS. 

Values for the vsim were taken from Rov er et alJ §002), 
computed from spectra collec ted a t the Haute-Provence 
Observatory (OHP) and bv lAerts et alJ J2004I) . derived from a 
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cross-correlation function analysis of spectrosc opic measure- 
ments with the CORALIE spectrograph. Mathi as et ail |2004) 
have published the results of a two-year high-resolution spec- 
troscopic campaign, monitoring 59 y Doradus candidates. In 
this campaign, more than 60% of the stars presented line profile 
variations which can be interpreted as due to pulsation. From 
this work, in which projected rotational velocities were derived 
for all 59 candidates, an additional value of vsim is obtained for 
HD 12901. 

The physical parameters used in the present work to delimit 
the error boxes for the star in the HR diagram are those given by 
TempLogG, derived from the Strtimgren-Crawford photometry 
and given in blold font in Table |5] Errors for these parameters 
will be assumed to be approximately 0.2 dex in log g and 200 K 
in T eff . 



Table 6. List of selected integer numbers sets associated to the 
observed frequency ratios of HD 12901. For each set, the cor- 
responding observed 3 bs is given in yuHz. An estimate of the 
spherical order I is given in column 4. 





hi 




"3 


{ 




h 


17 


27 


31 


1 


987.0 


h 


21 


33 


38 


1 


1202.4 


h 


21 


33 


38 


2 


694.2 


U 


26 


41 


47 


2 


860.0 


ts 


30 


47 


54 


2 


984.3 


k 


33 


52 


60 


2 


1087.9 



4.2. Modelling 

The evolutionary code CESAM dMorellll997l) has been used 
to compute stellar equilibrium models. The physics have been 
chosen as adequate for intermediate mass stars. Particularly, 
the opacity tables are taken from the OPAL package 
Jlglesias & Rogerslll99fJ) . complement ed at low temperatures 
(T < 10 3 K) by the tables provided bv lAlexander & Ferguson! 
( 1994). The convective transport is described by the classical 
ML theory, in which the free parameters cxml = lm/Hp =1.8 
and a mixed core overshooting param eter d„ v = l„„/ Hp = 0.2 
are considered (as prescribed by Sch aller et alJI 19921 for inter- 
mediate mass stars). The Hp corresponds to the local pressure 
scale-height. l m and d ov represent the mixing length and the in- 
ertial penetration distance of convective bulbs respectively. For 
the atmosphere reconstruction, Eddington's T(j) law (grey ap- 
proximation) is used. The transformation from heavy element 
abundances with respect to hydrogen [M/H] into concentration 
in mass Z assumes an enrichment ratio of AT/AZ = 2 and 
T pr = 0.235 and Z pr = as helium and heavy element pri- 
mordial concentrations. 

To cover observational errors in the determination of T e g , 
log g and [Fe/H], equilibrium models are computed within a 
mass range of 1.2— 1.6M Q and metallicities from -0.6 up to the 
solar value (Fig. @J. Theoretical oscillatio n spectra are com- 
puted with the oscil lation code FILOU (see lTran Minh & Leonl 
1995>. ISuareJ2 002). Eigenfrequencies are computed for mode 
degree values of up to t — 2, thus covering the range of ob- 
served g modes for y Doradus stars (up to g^). 

4.3. Locating HD 12901 in a J -logT eff diagram 
4.3.1 . Observed frequency ratios 

Following the steps described in Sec. [3] the ratios of observed 
frequencies are: 

A = 0.871, ^= 0.639, -^= 0.556. (12) 

On the other hand, all possible integer number ratios up to 
n — 60 are calculated. An error of 5 • 10~ 3 on ratios values is 
assumed in order to match the theoretical uncertainties found 



in Fig. [3] With these assumptions and considering only dipoles 
{( — 1) and quadrupoles (( = 2), as has usually been found in 
this type of star, all possible sets (~ 60) were reduced to six, 
given in Table|6] For each set, the spherical order t is estimated 
by comparing their radial orders with theoretical spectra for 
representative models (see Sec. l4.2l and Table. 1). Furthermore, 
it is possible to obtain the observed Brunt-Vaisala' integral 
(3 bs) through Eq. [2 constituting one of the main constraints 
of this method. 

4.3.2. The J-logT eff diagram 

For each valid set of integer numbers (f, in Table|6j there is only 
one associated 3 oos value, as deduced from Eq.[2 Selected U b s 
can thus be located on a theoretical U-logT e ff diagram. This is 
done in Fig.|5]for the selected f,-. Only t\, tj, and ts, are displayed 
since ti and t b give a Brunt-Vaisala integral value impossible 
to reproduce with the theoretical models (see Fig. [5}, and fg 
has a 1 bs impossible to differentiate from t\ in the 3-logT e ff 
diagram. The theoretical 3 is computed for a set of representa- 
tive models of HD 12901 (Fig.gJ in the range of 1 .2-1 .6M for 
three different metallicities ([Fe/H] = 0.0, -0.4 and -0.6). 

Models are computed within an error box that takes into 
consideration the typical accuracy for the calculation of physi- 
cal parameters. As discussed in Sec. 2, a numerical uncertainty 
in 3 of 1% is chosen. 

4.4. Toward the modal identification 

At this stage, the choice of good models is drastically reduced 
to those within the 3 b s -log Trff boxes obtained for each f, set 
given in Table |6] Nevertheless it is possible to provide further 
constraints on both equilibrium models and modal identifica- 
tion from additional observational information. 

Considering the observed metallicity for HD 12901 (see 
Table |3 and the corresponding error in its determination, only 
models with [Fe/H] =-0.4 and -0.6 are kept. This reduces the 
set of possible valid identifications to fi, £4 and ts with their cor- 
responding spherical order, I and 3 t, s values, shown in Table 6. 
This constrains the mass of models to the range of 1.2— 1.4M . 
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Evolution ^ 



3.90 3.85 3.80 3.75 

logT eff 



Fig. 4. Evolutionary tracks of models representative of HD 12901 in a g- T e g (in a logarithmic scale). The shaded area represents 
the observational error box as deduced from Table[5] From bottom to top, masses of models vary from 1 .2 to 1 .6M Q . Continuous 
lines represent evolutionary tracks of models computed with [Fe/H]=0.0 (the solar value). Dotted lines represent those computed 
with [Fe/H]=-0.4 and dashed lines are those computed with [Fe/H]=-0.6. 



Table 7. Charateristics of the selected best models obtained by comparing theoretical frequencies of deduced sets from represen- 
tative models with 3,/, ss 3 b s . Effective temperatures and luminosities are given on a logarithmic scale. Ages are given in Myr. 
For each selected model, the corresponding mode degree t is also given. 







t 


[Fe/H] 


M/M 


T e ff 


L/L 


log g 


x c 


Age 


P/Po 


A\ 




2 


-0.4 


1.2 


3.83 


0.52 


4.28 


0.50 


1990 


9.19 


Ai 




2 


-0.4 


1.3 


3.85 


0.72 


4.17 


0.40 


2100 


6.10 


A 3 




2 


-0.4 


1.4 


3.84 


0.90 


4.02 


0.26 


2090 


3.47 


Si, 


Ci 


1.2 


-0.6 


1.2 


3.83 


0.67 


4.12 


0.27 


3120 


5.36 


Bi. 


c 2 


1,2 


-0.6 


1.3 


3.83 


0.84 


3.98 


0.17 


2720 


3.20 


B 3 , 




1,2 


-0.6 


1.4 


3.83 


0.98 


3.88 


0.10 


2290 


2.20 



4.4.1. Theoretical frequencies 

Another constraint for selecting the theoretical models fitting 
the observations consists of comparing, for each set, the theo- 
retical frequencies corresponding to each «, with the observed 
frequencies. To do so, theoretical oscillation spectra have been 
computed in the range of g modes given by selected t\,U and 
?5 sets for models within the error (shaded) boxes of Fig. [5] 

In Fig.|6] the set of models verifing the integral value of t\ 
is depicted. All these models have thus a metallicity of -0.4. 
From top to bottom different masses are shown (1.2, 1.3 and 
1.4M ). In left panels the frequencies for the modes g26, #41 



and g47 (see chain t£) obtained with each model are represented. 
The horizontal lines are the observed frequencies of HD 12901. 
From these panels it becomes apparent that there is a model 
verifing fairly well the complete chain obtained by our method 
as well as the observed physical characteristics. The right pan- 
els show the relative error defined as follows: 



1 V bs 1 

The choice of the best model within a track can be carried out 
by searching for the model minimizing this relative error. This 
model is represented in Fig. |S]by filled circles. Repeating this 
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Fig. 5. J as a function of the effective temperature for representative models of HD 12901 with masses in the range of M — 1 .2- 
1.6M , for three different metallicities [Fe/H] =0.0, -0.4 and -0.6. Shaded areas represent uncertainty (3 b s ,log T e fr) boxes of 
selected sets t\, tj, and ts, as given in Table[6] 



procedure for the rest of the selected models in Fig.[3]it is possi- 
ble to restrict the complete sets of models to just nine. Seismic 
models with the same mass and metallicity but verifing modal 
identifications with different I are computed from the same 
equilibrium model. This reduces the total number of theoret- 
ical models fitting HD 12901 to 6 (those verifing the complete 
constraints concerning log T e ff, log g, metallicity, (n, /, <x) and 3 
for each observed frequency). In Table [7] the charateristics of 
the selected best models are displayed. 

This final set of models is depicted in the HR diagram in 
Fig. 9. As can be seen in Tabled this set is very heterogeneous, 
with very different hydrogen central concentration, age, log g 
and/or mean density. This is particularly useful when addi- 
tional constraints on models exist. In this context, information 
on the spherical degree of modes (from multicolour photometry 
(Mo va et alJ20 04) or spectroscopy) would be speci ally helpful 
in discriminating between models. In particular, in lAerts e t al. 

the three frequencies of HD 12901 are suggested to be 
I = 1. This would restrict the representative models of the star 
to B u B 2 and B 3 . 

Further constraints from the pulsational behaviour of 
y Doradus stars are likely to complete the method pre- 
sented here. Particularly, constraints on unstable modes con- 
sidering a convection-pulsation interaction can be obtained 
(work in preparation) fr om recen t theoretical de velopments of 
Griga hcene et alJ d2004l) and lDupret et alJ (12004 . 



5. Conclusions 

By using the first order asymptotic representation for the low- 
frequency g-mode eigenvalues of stellar pulsations, a method 
for estimating the radial and spherical orders as well as the 
Brunt-Vaisala integral is presented. The method is developed 
under the assumptions of adiabaticity and no rotation, and is 
based on information gathered from ratios of observed frequen- 
cies, similary to the large and small frequency difference for 
high-order p modes. It can be shown that the method becomes 
useful when applied to at least 3 frequencies. For each of these 
observed frequencies, some information about the I of each 
mode or the assumption that all of them have the same t is 
needed. 

The frequency method is specially adapted to 
y Doradus stars for two main reasons: 1) they oscillate 
with g modes in the asymptotic regime and 2) their internal 
structure (a radiative envelope between two convective zones, 
the inner core and other close to the photosphere) make it 
possible to use an asymptotic expression with only one scaling 
model-dependent parameter, the Brunt-Vaisala integral. 

The self-consistency of this method is verified by means 
of a theoretical exercise in which the accuracy of predictions 
given by the expresion obtained bv lTassoull ( 1198 0) for asymp- 
totic g modes is checked. These predictions are compared with 
theoretical frequencies computed for a model representative of 
a typical y Doradus star. 
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3.8325 



3.3320 



3.831 e 



3.8310 



fi t- « 



3.8325 



3.8320 



3. 83 1 5 



3.8310 



35 - 



i 20 



15 



[Fe/H] = -0.4, t 4 -1.3M e 



3.850 3.848 3.846 3.844 3.842 3.840 



[Fe/H] = -0.4, t 4 -1.4M Q 




3.8510 3.8500 3.8490 3.8480 3.8470 3.8460 3.8450 



1.0 -i 



[Fe/H] = -0.4, t 4 -1.3M e 



• - • • • 



3.850 3.848 3.846 3.844 3.842 
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2.0 - 



0.5 
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Fig. 6. Evolution of theoretical frequencies (left panels) (g26,gA\,gAi',( = 2) corresponding to the set t\ (see Tabled for selected 
1.2, 1.3 and 1.4M models with a metallicity of [Fe/H]=-0.4 and 3 t h = ^obsiU)- Filled circles represent such frequencies for each 
model within the photometric error box of Fig. 4. The observed frequencies fi, fn and /m are represented by continuous lines. 
Right panels show the error ej v of theoretical frequencies with respect to the observed ones for each panel on the left. Frequencies 
of best models are represented by black circles. 
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Fig. 7. Surface gravity as a function of the effective temperature of seleceted best models following the procedure given in Fig. [6] 
These selected models are represented by filled circles which are labeled in Table Evolutionary tracks of concerned models 
are also displayed: as in Fig|4] dotted lines represents tracks for models with [Fe/H]=-0.4 and dashed lines those for models with 
[Fe/H]=-0.6. The shaded area represents the observational error box deduced from Table[5] 



An application to the real star HD 12901 is also given. 
This star has been recently found to be a y Dora dus star and 
three oscillation frequencies have been identified (lAerts et alJ 
120041) . Nine possible mode identifications (n,() as well as the 
corresponding Brunt-Vaisala integral estimates are obtained. 
From standard constraints on models coming from fundamen- 
tal parameters (mainly log g, T e fr and metallicity), nine sets of 
representative models are selected. From these models, the- 
oretical frequencies computed for the nine mode identifica- 
tions obtained are compared with observations. This lowers 
the total number of representative models of HD 12901 to six. 
Furtherm ore, recent m ulticolor photometry results for this star 
given by lAerts et alJ J2004) 

suggest the spherical order to be 
i — 1 . This drastically reduces the number of valid models fit- 
ting the observations to only three. 

This method constitutes an important step toward the modal 
identification of y Doradus stars, especially in cases when there 
is no additional information, like that provided by multicolour 
photometry or high resolution espectroscopy. It is, therefore, 
a method well suited for the particular case of the white- 
light, very high precision photometry expected to be deliv- 
ered by the asteros eismology camera of the COROT mission 
(Baglin & Auvergn el 19971) . Nevertheless, improvements com- 
ing from the pulsational behaviour of y Doradus stars such as 
unstable mode predictions considering a convection-pulsation 



interaction, and additional spectroscopy or multicolor photom- 
etry information will complete the method here presented. 
Therefore, this method, together with the additional informa- 
tion mentioned above, represents a complete scheme for study- 
ing and interpreting the observational behaviour of y Doradus 
stars. 
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